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Abstract 

For the noncommutative torus T, in case of the N.C. parameter 9 = — and the area 
of T is an integer, we construct the basis of Hilbert space Tt n in terms of 9 functions of 
the positions Zi of n solitons. The loop wrapping around the torus generates the algebra 
A n - We show that A n is isomorphic to the Z n x Z n Heisenberg group on 9 functions. 
We find the explicit form for the local operators, which is the generators g of an elliptic 
su(n), and transforms covariantly by the global gauge transformation of the Wilson loop 
in A n - By acting on TL n we establish the isomorphism of A n and g. Then it is easy to 
give the projection operators corresponding to the solitons and the ABS construction for 
generating solitons. We embed this g into the L-matrix of the elliptic Gaudin and CM. 
models to give the dynamics. For 9 generic case, we introduce the crossing parameter 
r] related with 9 and the modulus of T. The dynamics of solitons is determined by 
the transfer matrix T of the elliptic quantum group A T , V , equivalently by the elliptic 
Ruijsenaars operators M. The eigenfunctions of T found by Bethe ansatz appears to be 
twisted by rj. 
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1 Introduction 

The development of soliton theory on noncommutative geometry is rather impressive 
(details of reference e.g. see review [|], g]). The paper || gives the solitons expressed by the 
projection operators. Witten Q gives the partial isometry operators. The paper Q uses the 
partial isometry operators to generate the solitons on noncommutative plane. 
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Recently, the solitons on the N.C. torus has attracted a lot of interest [||, [?], ||], ^, 10. [H], 12, 
13]. On the compactified torus [14], the duality, the Morita equivalence and the orbifolding 



have been studied [15]. 

The equivalence class of projection operator on torus is given by (l(| in terms of U\ and 
U2- The noncommutative algebra A generated by Ui (U\U2 = U?fJ\e. ) wrapping the torus 
has been given in [14, 17| in terms of the matrix difference operator. The matrix acts on 
a U(n) bundle V n with trivial connection while the difference acts on a U{\) bundle C. The 
Hilbert space 7i = V n ® C acted by A is given by the vector functions of real variables. The 
local covariant derivative operator v« acted on C is given also. Their commutator gives the 
constant curvature related to 9. Obviously the periodicity of wave functions on torus can 
not be given explicitly by real variable functions. The Schwarz space of smooth functions 
S(R) of rapid decrease real functions provides a bimodule between Aq and Ai under Morita 

equivalence [0, 12]. Starting from the Gaussian function e lny2 ^ e as the Schwarz function, Boca 



constructed the projection operators in terms of 9 functions. This projection operator 



satisfies the BPS like selfduality condition || 18], thus gives constant curvature. But Boca 
obtained the explicit expressions only for the case with modulus equals 4, moreover, it is in 
terms of the products of two 9 functions depends seperately on U\ and U2, so the symplectic 
structure of noncommutative torus is unclear. Gopakumar et al [1C] starting from the same 
Gaussian function but orbifolding the double periodic multisoliton solution on the N.C. R 2 
into a single soliton solution on the torus. They employed the kq representation [|H], 2C] 
on a dual lattices for conjugate variables k and q, which provides a basis of simultaneous 
eigenstates of commutative U\ and U2 with 9 = 2irA, where A is an integer. Thus they 
succeeded in constructing the soliton on the so-called integral torus with double periodic wave 
functions. The noncommutative symplectic complex structure on torus appears explicitly and 

the corresponding Weyl-Moyal transformation is realized as double series of and U 2 A ■ But 
since the loop U\ and U2 become commutative in this degenerated 9 = 2nA(e ld = 1) case, 
they obtained just a unique projection operator, there is neither the covariant derivative 
corresponding to the local translation nor the ABS like partial isometry operator. Moreover, 
it is not obvious how to generalize into the case of generic 9. 

This paper is just atempted at first to give explicitly more projection operators, the local 
derivative operators E a and the soliton generating ABS operators by orbifolding such integral 
torus T into a T n ; secondly to realize the noncommutative algebra Ag in case of generic 9 
explicitly on a Hilbert space of meromorphic functions. We try to relate various algebra on 
noncommutative torus, the Z n x Z n Heisenberg algebra with the elliptic su n (T) algebra and 
the elliptic Sklyanin algebra with the elliptic quantum group, then furthermore consider the 
dynamics of the noncommutative solitons by relate it to the well known integrable models 
such as the elliptic Gaudin models with the ellipitc Calogera-Moser models and the elliptic 
IRF (interaction aroud face) models with the elliptic Ruijsenaars-Schneider models. 

In the next section, after review the construction by Gopakumar et. al [1C] for the soliton 
on the integral torus T, we show by n x n times orbifolding this T into T n , that the one 
dimensional trivial Heisenberg group i.e. the double periodicity under wrapping by Ui on T 
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in paper |10| will be refined into a Z n x Z n G-^(n) generated by loops Wj(= [/>) around % 

In section 3, we use the 9 functions of the central "positions" z% of multisolitons to give 
an explicit expression for the basis of the Hilbert space 7i n . The A n algebra turns to be 
the Z n x Z n Heisenberg group Gy_(n) acting on 9 functions |23[ |. The local infinitesimal 
translation operator E a is realized as the derivative operator generating the su n (T) algebra 
on torus. The E a behaves covariantly under gauge transformation by Gy_{n). We establish 
the isomorphism of the algebra A n and the su n (T) by applying them on TC n , they both 
become identical to the G-^(n). Then it is easy to construct the ABS and the projection 
operators corresponding to multiple solitons. 

Later, in section 4, we embed this su n (T) derivative operator as the "quantum" operators 
in the representation of the transfer matrix (Lax operator) of the elliptic Gaudin model i.e. 
L(u) is obtaned from the classical Yang-Baxter r-matrix r(u) = J2a w a(u)I a <8> la acted on 
V n ®V n space by replacing one V n into a "quantum" space, i.e. replacing the Heisenberg 
matrix I a on it by the su n (T) operator E a . The r-matrix is nondynamic. It is depends 
on the spectral (evaluation) parameter (u) through the meromorphic sections w a (u). The 
double Heisenberg properties if w a (u) and hence of r(u) ensures L(u) to be a section on a 
twisted su(n) bundle. Then we gauge transform the Gaudin Lax into that of the elliptic 
Calogero-Moser models. This Lax is equivalent to the elliptic Dunkle operators representing 
Weyl reflections. The trace of the quadratic of L gives the Hamiltonian of the C. M. which 
is assumed to give the dynamics of the N.C. solitons on the brane. As in section 3 we have 
shown that the isomorphisim of N.C. Z{ and Z{ to the Zi and local translation di on H in N.C. 
1Z 2 case, under orbifold 1Z 2 /Z (g) Z becomes the isomorphism of the Z n x Z n Heisenberg W a 
in A n to the E a of su n (T). So the correspondence of noncommutative <9j(~ z{) and Z{ in E a 
to dynamical pi and qi in L(u), natrually endow the N.C. torus with symplectic structure. 
Thus naturally the solitons on branes satisfy the dynamics of Calogero-Moser systems. 

Then we turn to the case of generic 9. Here, at first in section 5, we find that, the 
difference representations of the Sklyanin algebra, where the difference operation is to shift 
the solitons position Z{ by rj, will give a constant Heisenberg matrix on the functions in V n 
introduced in section 3. The chracteristic determinant of this operator gives the Ruijsenaars 
operators. Since the Ruijsenaars- Schneider model in the limit r\ — > approaches the CM. 
model, it is reasonable to assume that the dynamics of the soliton on 9 generic torus is given 
by using these Ruijsenaars operators as the Hamiltonian. 

But this Sklyanin algebra on the vector space V n gives only the constant matrix part 
of the trivial su{n) bundle. It is not sufficient to show the whole su{n) (8) u{l) operations 



as [14, 15]. So in section 6 we realize the full algebra as the elliptic quantum group E TjT) 
constructed by Felder. It is well known that the E T „ is equivalent to a S Tt „, their evaluation 
representations is related to the elliptic IRF models. At first we give a factorized form of this 
representation, which may be gauge transfromed into the Sklyanin algebra and further we 
relate its transfer matrix to the Ruijsenaar operators M. These representations of S Tt „ are 
the to be found matrix difference operators which realize the corresponding operators given 
by [14, 15], but instead of acting on real vector functions, now it acts on meromorphic vector 
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sections. Then as Felder and Varchenko |21|], using the matrix difference operator M to 
represent the transfer matrix, we apply the Bethe ansatz to find the common eigenfunctions 
which actually twisted by r\ for N.C. f/j and realizes the Hilbert space |L4|, 15] in the form of 
the vector meromorphic sections (sheafs). 

In the last section, we shortly describe the subjects which will be investigated later. 

2 Solitons on the "integral torus" |TD|| and its further orb- 



ifolding 



In this section we at first shortly review the result of paper [10] for solitons on nocommu- 
tative space R 2 : [£1,^2] = iQ. This R 2 has been orbifolded to a torus T with periodicities 
L and tL. 

The generators of the fundamental group of T are: 

C/l = e~ 1 ^ 1 , U 2 = e^ 73 * 1 - 71 * 2 ), (1) 

where I is the normalized length I = Ly/Q, y l = (* = 1 5 2) and 

U X U 2 = U 2 U ie - iT212 (2) 
As in fpLOTl , let us consider the case T #^- £ N (or Z+) i.e. the normalized area A = = B 



the flux passing through the torus is an integer. The projection operators constructed by [1C] 
has its image spanned by the lattice of coherent states 

C/f^0),(ii,j 2 eZ 2 ) (3) 

here |0) is the Fock Bargmann vacuum vector: a|0) = 0, where a = -^{y 1 + *y 2 )> °^ = 
-^(y 1 — iy 2 ), [a, at] = 1. They found a particular linear combination: 

\^) = J2c jl , j2 Utut\0) (4) 

h ,32 

that satisfies 

m^ui 2 m = 8 n0 5, 2 o. (5) 

Then the projection operator becomes 

31,32 

P 2 = P. (6) 

Since dH), we have UiPU^ 1 = P, (i = 1, 2), i.e. P actually acts on T = with periods f/j. 
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To find the solution satisfying the orthonormality condition (0), the wave function of 
has been found [10] in terms of the kq representation [19, poj ]: 



\kq) 



' J_ P ~iri(y 2 ) 2 /2r 2 



(7) 



where |g) is an eigenstate of y . y \q) = q\q), such that 



JIT2V 1 



\q) 



Jlqr 2 



\q), 



\q) = \q + l). 



Obviously, in case A being an integer, \kq) is the common eigenstate of U\ and E/jj: 



Utlkq) 



-ilk 



\kq), 



U 2 \kq) 



Jlr 2 q 



\kq). 



(8) 



(9) 



Remarks: From the completeness of the coherent states [E2|, we know that \kq) is defined 
in the double periods: < k < < q < I, and constitute an orthonormal and complete 

basis for the Hilbert space TLa = Tt^/A, if A > 2. In case of A = 1, if we choose just one 
coherent state for any unit cell, then all the state excluding the vacuum state maybe taken 
as a complete and minimal system of state. Bars [ 13 [ recently has discussed the solutions 
with size smaller than this extremal one. 



By using \q) = exp(^^ — l~ ^1 + we find the wave function 

(Mo) 



tV 2 



(10) 



where r = T\ + it 2 - Then by modular transformation of ( [TQ| ) or directly from (|7j) by using 
the Poisson resummation formula we obtain as in [10]: 



C o (k,q) = {kq\0) 



— r^=exp(- — k +ikq)6 m ( , - . 

Kiy/l 2it 2 I A 



(11) 



The orthonormality condition (||) becomes 



^1,0^2,0 



dk I dqe 



^ rl 



-ijllk+ij 2 lr 2 q\ 



9(k,q)\' 



A-l 



dk / dqe 



-ijllk+ij 2 W 2 q\ 



%k,q)\ 2 Y: 



n=0 



C (k,q + n 



(12) 



Thus, paper |TJJ find that 

(kq\iP) =V(k,q) = c{k,q)C (k,q) 



C (k,q) 



27r/AEn=o\Co(k,q + nl/A)\< 



(13) 
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Now we turn to the further orbifolding. Let W\ = U" , then 



i , rk' 2 i„„ ,q + ^ t 



d(k,q) = (fcg|Wi|0) = (k,q + -|0> = _-exp(- — + ik(q + -))0 O i (— r^, -r), (14) 



n 



7T-i 



here we have chosen A and n relatively prime. 
Similarly 

Cp(k,q) = (kq\W^\0) = -i^exp(-^+*%+^))V(^^ + ->T) 

' : exp(-^ + ik(q+^l))e oE (z,^), (15) 



here z = yi + %2) the map k,q to z is given |10| by the Weyl Moyal transformation. The 
effect of W\ is to shift z to z + ^ , or equivalently shift the characteristic /? of 6 a ^ to /3 + - . 
Constructing 

= ,f^' g> , ■ (16) 

then l^+n) = e -4 ^!^^) and \^ a ) (a = 0, 1, • • • , n — 1) are linearly independent. 
Let 

Pp=Y. Utui^pWpp-^U^, (17) 

31,32 

then 

Pp = Pp (18) 
W^Pp 2 W^=Pp l+ p 2 . (19) 



As in paper [15|, We may also orbifold the torus along the direction n times. Let 
i — 
W 2 = U 2 n , then 

C' a (k,q) = (kq\WZ\0)=e i ^ q (k+—cx,q\0) 

n 

-j t 2 (7-1- — 

= eM-^+i(k + -lr 2 )q)9^ A0 C— ^1). (20) 

Obviously, W2 shifts z to z + ^, i.e. shift the a of a)( g. 

Subsequently, we may construct \tp' a ) and P' a as (16) and (|i"7|) . But the set P^, (a = 
1, 2, • • • , n) are not independent from the set Pp, (/3 = 1,2,---, n), and in either basis the W\ 
or(and) W 2 matrix are not constant. Actually, the target space Hf of P on total torus T 
has been subdivided into 7i n on torus T n as described by [15|. One may find a basis of 7~in 
such that W± and W 2 become the n x n irreducible matrix representation of the Heisenberg 
group 

W1W2 = W 2 Wiuj, oj n = 1, W? = W£ = 1. (21) 
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In next section we will construct explicitly this basis in terms of Q a p functions. But 
before that, let us compare to the case of N.C. plane, that the generic n solitons solution 
is Iir=i( a ^ — z i)|0), with n soliton centers Zj. Similarly, we will introduce the z$ for the 
location of the centers of n-solitons solutions on the torus T, such that the moduli space of 
the n-soliton solution is T® n / S n . 



3 Noncommutative algebra A n , Hilbert space H. n , Heisenberg 
group Z n x Z n and su(n) algebra, in case [£7i, U2] = 

3.1 The Hilbert space 



As in [15], usually the Hilbert space Tij- can be written as the direct product of a su(n) 
trivial bundle V n and an U(l) line bundle C: TLr = V n (& £. On the V n acts the W\ and 
W2 matrices, on the C acts the covariant defference operators Vj (the notation of |14|), and 
Wi ® Vj = f/j. Thus, the operators U( are matrix difference operators acting on vector 
functions v a , (a = 1,2, • • • , n). But now it happens that in case of commutative U\ and U2, 
this V n is identical to the whole TL n on subtorus T n in the following way. 

The basis vectors of the Hilbert space H n are 



£F_ a>6 ,(a = l,2,--., 



b=i 



TP = 7? - p«7T««2 



n^i 



1 n 



fc=i 



here a = (ai, 02) G x Z n , and 



a a (z) 



1 I Oil 

2 "1" n 



(In the following, we use the modulus r to represent the J in section 2) 



(22) 



Remark: The ^ (13) has been normalized, but it is not entire, i.e. it has poles and is 
nonanalytic. Our basis v a is entire, they span the space of weight n quasiperiodic functions 
p3| , i.e. the space of n soliton sections of some su(n) bundle. 

Now we turn to show that the W" 1 
generates the algebra A n on T n . 



W2 2 with (0:1,02) € Z n x Z n acting on this H n 



3.2 Noncommunitative algebra A n on fuzzy torus T as the Heisenberg 
Weyl group Z n x Z n 

From section 2 we know that the effects of the noncommutative Wilson loop [24, |25| W\ = 

Ui and W2 = (U 2 n ) on the i-th soliton is to translate its position, from Zi to (zj+-r— 5i n r) and 
(zi + - — 5i n ) respectively, or equivalently shift all Zi by ^ or ^ mod the torus T, furthermore 
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equivalent to shift the coordinate origin u by ^ of ^ in oposite direction. Substituting in (22), 
we find 



rt-i 

^W\T-.( n ) 



= C[[T^)Tr T v a (z 1 ,---,z n ) (23) 



i=l 



Vo-l(^l, • • • , ^n) = Wiv a (zi, ■■■ , z n ), (24) 

«o(^l + -,••• H 1) (25) 

n n 

(-l^e^vaizx, ■■■,z n )= W 2 v a ( Zl , • • • , z n ), (26) 



where 



= e™ ^^/(zi, • • • , zi + or, • • • , z n ). (27) 
Using basis v a (p2[), we have transform the action of the N.C. Wilson loop of T n from the 



shift in functional space form (|23| ) and (25) into the matrix operator form (f24|) and (^3), i.e. 

(Vi)«6 = S a+1}b ; (V 2 )ab = S ab uj a , u = e^=e i6n (28) 
As expected this Vi and V 2 satisfy the relations 

vtVi=Vjv 2 = l, V? = V2 = 1, V 1 V 2 = V 2 V 1 e 27Tie \ (29) 

here 6 n = ^ is the noncommutative parameter for the n 2-branes [15]. 

Remark: There are some subtle and crucial points in our paper which usually cause 
confusion. It seem worth to be stressed and clarified here. Originally for Tit on total T we 
have the commuting U\ and U 2 with 9 = 2irZ. After orbifolding, it is clear [14, 15 1 that 
T~Ct = V n ® C is acted by Ui = Wj <S> V^" 1 , where the matrices Wj satisfy W\W 2 = W 2 W\uj, 
while the difference operators 

Vr 1 ^ 1 = V^V^e -™ = il 2 l Wl l e~ 2me -. (30) 

Generally is represented by Z n x Z n matrix on trivial su(n) bundle V n , while V" 1 is the 
translation of the origin of T: u = — ► u = — 21+221 ? correspondingly, to shift the center 
of mass of n solitons z c = ~ Ya=i z i by ai +" 2T , If we restrict to consider the H n = V n on 
subtorus T n , then we may either realize the W-* in operator form i.e. in terms of Z n x Z n 
matrix ( p6|) or realize it in functional formalism i.e. in terms of difference operators ( |23|) 
(|25|). These two forms is related by the Weyl Moyal transformation on torus. But it happens 
that the V a (f c — ► z c + 2i±p22:^ yj e u a ma trix transform on our basis function, so it is 
equivalent to the operator form of W a . Thus for the total torus the W a and V" 1 cancels, 
U a = W a <g> V" 1 are commutative, as the (6.54) in [15]. 

The functional and the operator forms are related by Weyl Moyal transformations on 
torus (see Appendix B of [0]), functions v a behaves as the coherent state functions on torus. 
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On the space H n = {v a \a = 1, 2, • • • , n}, the noncommutative algebra A n generated by W\ 
and W2 has been truncated and becomes a n x n dimensional unital C* algebra with n 2 basis 
W a and the define relations 



W^W? 2 



W a/3 w a^ = ^ w aP w a'f3' = w a' (3 1 wa /3 -(3a 1 = yya+a' ,f3+/3' ^-fia' ^ (W<*0)n 



1. 



(31) 

Here and hereafter to represent in matrix (operator) form the basis W ai,OL2 of A n we introduce 
the usual Z n x Z n matrix (I a )ah = (4i,a 2 )a6 = 5 a+ait a 2 u> ba2 |23|]. This algebra in difference 
operator form is the Heisenberg group G?{(n) Z n x Z n in ordinary 9 function theory 23 j, 
including both the shift + ^ and the change of phases. 



3.3 su(n) algebra g on torus 



It is shown [26] that the level / representation of the Lie algebra sl n (T) on the elliptic 
curve T can be written as following: 



£ Q = (-ir^o) £n 



and 



3 k^j 



En 



l_ a 'a( z ji) _ q 



E^ 



(32) 



(33) 



here a = (01,02) £ Z n x Z n and a ^ (0,0) = (n,n), z 



■jk 



Zk: 9 = & Eq commutes 



<) 



with E a , sl n (T) includes only the E a with a^O. After a complicate calculation, we obtain 
the commutation relation: 



■c*27i 



-ai 7 2) E 



a+7) 



or in usual sl(n) with i,j label the Chan Paton indices basis, let 

Eij = ^ {I )ijE a , 



then 



[Ejk, Ei m ] — Ej m 5ki — E[ k 5 



jm- 



(34) 



(35) 



(36) 



Remark: The representation (22) and the commuation rules ( |34| ) can also be obtained |2 



from a quasiclassical limit from the representation of the Z n x Z n Sklyanin albebra 27]. 



Sklyanin and Takebe |2£] give the elliptic sl(2) by using double periodic Weierstrass functions. 
The high spin I representations is given also. In this paper, it is restricted to the I = 1 
representation of sl n (T) by holomorphic sections on T® n /S n . 
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3.3.1 Automorphism of En G SU2{T) by noncommutative gauge transformation 

W a G A 

Since the Wilson loops W\ and W2 acting on the noncommutative covering torus T is to 
shift Zi to (zi + ^ — (5j n T) and (zj + ^ — <5j n ) respectively, so E a in (|22"|) will be changed into 



W 1 -(E a (z i )) = W 1 E a (z l )Wi L = E a (z l + --5 in T)=u J - a2 E a (z l ), 

n 

W 2 -{E a {zi)) = W 2 E a (z l )W^ 1 = E a {z i + --8 in )=u)^E a {z i ). (37) 

n 

or more generally 

W frjh . Ea = w^E^W^ 2 )' 1 = uj ail32 ~ a2Pl E a . (38) 



This could be compared with the matrix model on noncommutative torus [14], where 

UiXjUr 1 = X, ■ <y,f>-Rj 
or more exactly with the covariant derivatives 

Ui Vi Uf 1 = <% V j 

i. e. 

E a exp(- V l) ai exp(- V2 r. 
n n 

3.3.2 Isomorphism of su n (T) and A n on 7i n 
let E a £ g to act on v a , we find that 

#a^a = ^2(Ia)baVb- (39) 

As in section 3.2 we already learn that W a v a = J2i^a)ba v bi so ° n W n we establish the 
isomorphism 

su n (T) =>- .4 

£ Q — ► W a (40) 

Obviously this is correspondent with the noncommutative plane case, where one have 

iejid Xi F(x) => [x,*f(x)\ (41) 

and acting on Fock space 

8 z ^a, 8- z ^(J (42) 

Compare eq.(pl) with eq.( pq) and eq.(42) with (|3~2"1) and (|39|), it is easy to see that the 
infinitesimal translations on plane R 2 corresponds to su n (T) and the algebra ^4r,2 correponds 
to the algebra A n . Meanwhile the adj operation for _4r2 (fil] ) changes to the Adj operation 
(gauge transformation p3, |25|) for A n (ET~ 
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3.3.3 Vi and E a as generators of the Weyl reflection group 



Kac [29] has shown that the 9 function with characteristic a = («i,a2) £ Z n ® Z., 
transforms under affine Wely reflections as the Heisenberg group generated by Vi (|23|), (|25|). 
As we have shown in (24), ( p6| ) and ([M]), both and E a act on w a as I a , operators V\ and i<a 
as the coxeter element I\, operators Vq, and E<i as the cyclic element I2; operator / j) 
permutes i,j i.e. gives the reflection fij reversing the root ej — Cj. 

Here, let us stress the relations of N.C. algebras TC n , A n and su(n) with the solitons. 
Obviously, since in (|22|) each term F a are symmetrical with respect to permutation S n of 
the centers Zi of solitons, so the vectors v a is defined on the module T® n / S n . The operator 
Vi shifts each soliton Zi from the i-th covering of T over T n i.e. from the i-ih brane to the 
next one, or equivalently changes each brane to the previous one. The operator V2 shifts the 
U{l)® n phase of i-th brane by uj % . The Chan Paton su(n) of n branes is generated by Eij. 
Acted on 7i. n , both Vj and E^ generate the Weyl reflection. The offdiagonal elements give 
permutations of the branes, the diagonal elements generates the phase shift of soliton wave 
functions. 

3.3.4 Projection operators and soliton generating ABS operators 

In the operator formalism, the projection operators becomes 

lE^( W ^P,^, ,43) 

and the ABS operators is simply 

E w ^ Wl = Y: |Fa+l)(K| r- (44) 

a ((V a+1 \V a+ i)(V a \V a ))2 

The explicit functional expresions may be obtained by Weyl Moyal transformation adopted 
to the torus as described in the Appendix B of [TLT 



4 Affine algebra on N.C. torus, integrable elliptic Gaudin 
and Calogero Moser models [|2B| 



Now we will show that the cotangent bundle for N.C. elliptic su n (T) ( |32"|) with twisted loop 
Wi realizes the elliptic Gaudin model on N.C. torus, which is related to critical level twisted 
su(n) WZW i.e. A n _\ affine algebra on N.C. torus. Then use the elliptic Vandermonde 
determinant to gauge transform it into elliptic CM. model which is the Hamiltonian reduction 
of the cotangent bundle for the algebra of the semidirect product of the Cartan torus and 
the Weyl group |3(J . 
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4.1 N.C. Elliptic Gaudin Model 



The elliptic Gaudin model on commutative space [28] is denned by the transfer matrix 
(quantum lax operator): 

Uj{u) = w a(u)E a (I a ) i: j, (45) 

a^(0,0) 

where 

e'(0)a a (u) 

Wa{u) = tttt 7-r, (46) 

a a {0)a {u) 

E a and I a are the generators of su(n) (or A n -\ Weyl) and G-n(n) respectively. Using the 
general defining relations of su(n) (^) and Gj{(n) pi]), we find that the commutators of L 
(quantum version of fundamental Poisson bracket): 

]L (1) (ui),L (2) (u 2 )] = [r( 1 ' 2 )(n 1 ,n 2 ),L«(n 1 )eL( 2 )(n 2 )], (47) 

where the classical Yang-Baxter matrix 



^ «>«(«)(/«)? ® (OS (48) 

a^(0,0) 



is antisymmetrical 

*#(«) = »#(-«) (49) 
and satisfies the classical Yang-Baxter equation: 

[^12(^12), ri 3 («i 3 )] + [ri 2 (ui 2 ), r 23 (« 23 ] + [^13(^13), r 23 (u 23 ] = 0, (50) 

where Thus this system is integrable. Indeed, we can find from the quantum 



determinant of L the complete set of Hamiltonian. Some examples may be found in [28, 31 



32]. It was obtained also as the nonrelativistic limit of the Ruijsenaars-Macdonald operators 
which will be described in next section. As a lattice model, the common eigenfunctions and 
eigenvalues of Gaudin model is solved in terms of Bethe ansatz [33], which has been expressed 
by the conformal blocks of the twisted WZW models on the torus [34, p5| . 

To relate these well known results about this usual Gaudin model with that on the fuzzy 
torus T n (section 3), let us substitute the differencial representation (p2]) of su n {T) E a into 



the Gaudin L ((45|), then it turns to be in a factorized form [26]: 
Eo+ w a (u)E a {I a )) = Liu)) = Y,Hu,z) i k (l)- 1 (u,z)jdu - ij^^i^^k^i^^j, 

0^(0,0) k k 

(51) 

where the factors are the vertex face intertwiner 

71 n-1 

4>{u,zfj = 9i_i i(u + nzj - 2J z k H —,nr). (52) 

k 1 
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We can also show that ( jsT] ) is a representation of L by subsituting it directly into (47). Now 
the expression of r becomes 




where 6 l (u) = 9i_± i(w,nr). It is easy to find that it is Z n (g) Z n symmetrical r(u) = 

2 ™ ' 2 

(J Q x I a )r(n)(/~ 1 tg) I" 1 ) and turns to be equal to (|4g|), i.e. actually the factorized L* (|5l|) 
represent fl4"7|). 

The intertwinner " factor" s ( |52|) , now intertwines the Chan Paton su{n) index i for the 
brane (vertex model index) with the dynamical indices (face model indices) of dynamical 
soliton position Zj on the world sheet. 

Feigin et. al [33] has established the relation between the critical level su(n) WZW models 
and the rational Gaudin model, that is the L(u) and R(u) as ([45|) and (f48D, but instead of 
(|46|) now with a rational w^(u) = -, C. Kuroki and Takebe |32| find the same relation 
for the elliptic case, but with the WZW twisted, i.e. defined on a twisted bundle g tw : 

g tw := (C x g)l ~ 

where the equivalence relation ~ are (u, g) ~ (u + 1, hgli 1 ) ~ (u + r, ^gl^ 1 )- Correspond- 
ingly the w a (u) ( [46|) is a meromorphic section of the twisted bundle with single pole at n = 
and residue 1. While wrapping around different one cycle on the base torus the global gauge 
transformation for sections in the fibre are noncommute, but the base torus itself has the 
complex conformal structures only, it has neither metric nor symplectic structure, thus re- 
mains to be commutative. For our N.C. torus we have both z = y\ + iyi and z = y\ — iyi 
which are N.C, it causes both the global N.C. W\,W2 and the local N.C. z,z. 

For the usual Gaudin [3l]] on a marked commutative torus, the Zi is the marked points, 
i.e. the poles in L = J2i J2 a w a {u — Zi)E a I a . The pi is the conjugate momentum: classically 
{Pi, Qj} = Sij, quantumly [pj, qj] = 5{j. So, the L-matrix is endowed with dynamics, by finding 
Poisson (equivalently Konstant-Kirrilor) brackets at first, then by quantized the pi ~ 
to find the matrix differential form of L ||26||. While now for the Gaudin on noncommutative 



torus, the z\ are the center (position) of the z-th soliton, <9j as its infinitesimal translation is 
equivalently to the Here we should stress that the key point to translate the usual inte- 

grable models on commutative space to that of N.C. space is to notice that since in N.C. space 
[z~i,Zj] = 5ij, so N.C. plane automatically become the symplectic manifold corresponding to 
the quantum phase space, pi — > Zi, qi — > Z{ i.e in the holomorphic Fock Bargmann formalism 
[z,f] automatically becomes d z f just like p — > J^. Apparently the complex structure with 
given metric gives the symplectic structure, moreover it is automatically "quantized", as the 
geometrical quantization. 

Let us show that the dynamical equation induced by N.C. this way really gives the inter- 
action between N.C. solitons pi, by studying the potential of CM. model. 
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4.2 Elliptic Calogero Moser model and its equivalence with elliptic Gaudin 
model 

The elliptic CM. model is defined by the Himiltonian: 

n 

■i=i 

where p(z) = d 2 a(z). This quadratic and other higher Hamiltonian are generated by the 
Krichver Lax matrix 

L K r(u)) =8,5) + {I- fyJg a< )\ + f Zji \ . (55) 
J a(u)a{Zji) 

By the Poisson transformation (classically d{ ~ pi, Z{ ~ qi) 

Pi — >Pi - o— ln n *(q), (56) 
oqi 



here 



n(g) = n Va = --^(o), (57) 

.jr n 



and Lxr('u)*- becomes 



Lcm(u)} = ( Pi - L^-lnU(q))S) - -a'(0)(l - # ^V?' V (58) 

This may be further gauge transformed into the factorized L (^) of Gaudin model by gauge 
transformation matrix 

0(uig ,} = ? ^L. (59) 

It is well known that CM. model gives the dynamics of a long distance interaction between 
n-bodies located at Zi, (i = 1, 2, • • • , n). Now we will show that it is the interaction between 
n-solitons on the fuzzy torus while Zi becomes the positions of the centers of each soliton. 
According to paper |n|, for N.C multisolitons the potential term is argumented to be the 
Laplacian of a Kahler potential K, which is the logrithm of a Vandermonde determinant. 
Actually we have 

Ep(*;j) = E # kg n*(*i-*jO = E ( 6 °) 

i^j i j^k i 

= det(cj){) = a(nu+ - > )l[a(z l3 ). (61) 
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The variable u of the marked torus is the so called spectral parameter or evaluation parameter 
of Lax matrix Lxr(w)*-. 

As [|36f| , replacing the "spin" index i,j elements Eij (|35| 



in the Lax matrix Lxr (u)j of the 
permutation Sjj in A n -i Weyl reflection, LxrW; 



CM. model by soliton exchange zi «- 
becomes like the Dunkle operators [37]. This "spin" index exchange equivalence with particle 
(soliton) exchange is obtained if restricted to be acted on the space of wave functions which 
is totally symmetric in both spins and positions. This is the A n -\ Weyl symmetry for the 
su n (T) Chan Paton index of the branes and for the positions z%, Zj of solitons. 

Here the determinant of the vertex- face interaction is the determinant of our gauge trans- 
formation matrix from Gaudin to CM. We should mention that its u independent factor 
Y[(z) is the Weyl antiinvariant ground state wave function of C. M. It is the "phase func- 
tions" part of the conform block of the twisted WZW model on elliptic curve related to the 
Gaudin model |03|, |34|, [3q ]. Dividing by this phase the KZB equation of the elliptic Gaudin 

to the heat equation associated to the elliptic C. M. equation. 



Model reduces [34, 35 



5 The Z n x Z n Heisenberg group realized by the Sklyanin al- 
gebra S TtTI in case of general 9 

In previous section, we have show the utility of ellipitc Gaudin-Calogero model for the 
N.C field, brane and QH fluid in case of rational 8, it is reasonable to further find some 
kind of N.C. solitons for the generic, i.e. deformed 8. Acturally, there is a lot of hints, 
that the Ruijsenaars-Schneider (RS) model, as the deformed version of CM., descrbes some 
brane configurations. The elliptic RS and the infinitely related IRF etc. may be the hopeful 
candinate. 

Firstly let us remind the physical meaning of unrational 8. If we fix the area of the base 
torus, but change the field B ~ , such that the passed B flux quantum ^ (" normalized 
area") generally is not an integer, then the commutator XJ^UylJ^U^ 2 = e l ® becomes generic, 
unrestricted to the root of unity. This implies also that the ratio of the area for the Narain 
dauls SL(2Z) (g) SL{2Z) (modular for base T, Morita for target) is generic, i.e. unrestricted 
to be rational. We conject that it is related to the characteristical parameter rj of IRF model 
etc. i.e. we assume that 8 = rjr mod modular (Morita) transformation. 

Now instead of simply shift by - and ^ under the loop action Wi (^) and ( p5|) or the 



su(n) covariant differencial operator E a , we introduce the difference operator S a p7| , 38 1 as 
following 

S a F = ^2S^F(z 1 ,z 2 ,---,Zj-r],Zj + i,---,u n ), (62) 
j 

where 

*W = e^a a (v) n n ? ( ?7 nZfc + ^ (63) 
jg. ao(rvn)ao{nzj - nz k ) 
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Remark: The commutative relation of define S a is given in [B7]] for n = 2 and j38| for general 
re. The difference representation of S a is constructed and proved via Belevin's i?-matrix by 
1 27, and then proved directly by [||] and |39| ] 
Fundamental representation 

Applying S a with level / = 1 on V a G W n (p2|), we get 

(64) 
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Actually, acting on 7i n , the fundamental representation of S a £ S T)TI becomes the Heisenberg 
Z n x Z n matrix I a . We prove this by three lemmas. 
Lemma 1. The SqFq is independent of 77, 



SqFq — Fq. 



(65) 



From pi), (pi) and (E9) we get 



ao(nr))S Fo = ncr (r?)^ 



n 



cr (nzj - re2 fc + 77) 
cr (nz,- - nz k ) 



■cr (nzj - ^2 z k - (n - 1)77) J| - z fc + 77). 



(66) 



As a function of r/ with given Zi we find that the right hand side (r.h.s) of ( |6"6| ) has n 2 zeros 
at 77 = \[aiT + a 2 ). 

In the quasiperiodic phase of the r.h.s under shift 77 — ► 77 + r, the 77 dependent part is 
exp(— 27rm 2 77), that is the r.h.s has exactly re 2 zeros as the function do (77,77). Thus SqFo is 
independent of 77. Take 77 — ^ limit we find (|65|). 
Lemma 2. 

5 00 F 7 = F T (67) 

This can be obtained from lemma 1 by further applying Heisenberg shift z\ = z\ — ^(7^+72) 
on lemma 1 i.e. on both side of (|66|) as following: 



oo(wzj - nzfc + 77) . 
t7 (n;z,- -nz k ) ■ 



k+j 

a^(nzj - ^2 z k ~ (n - l)ry) J| cr 7 (re^ - ^ z k + 77) 
-o"o(n7/) Y[( nz i -y2 z k), 



(68) 



or simply 



(69) 
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In fact, it is independent of which Z{ only that is shift to Z{ — h(ji T + 72)- Since in the 
original and shifted F a , the nzj — J2 z k are weight vectors of A n -±, all will be shifted by 
(71,72) in the coset (weight lattice)/ (root lattice). While in S^o the root n(zj — Zk) with Zi 
(without Zi) changes (unchanges) by integer (71 r + 72). So there appears only a phase shift 
for each do, which cancels in total. It is easy to see that this is equivalent to the Heisenberg 
transformations (^), (^) and (|37]) on TL n where we have shift all the zi at the same time to 
the next coverings etc. 
Lemma 3. Transform lemma 2 by 

rj — > rj + — (air + a 2 ) (70) 
n 



we have 



q TP 27ria2(7i+ai)/n p 

S a F p = cu a ^Fp_ a (71) 



Really substituting (|70|) into (|67|), we get 

a a (nzj — nzk + f]) . 



0- o («^j - nz fc ) 



a^ +a (nzj - ^2 z k - (n - l)rj) J| a 1+Q (nzi - ^z k + rj) 
^a {nr,)e^^ +a ^ \{ a^nz, - £ z k ). (72) 



For higher representation, the fussion rule I = l\ + I2 for S a has been given in [4C]. The 
representation space basis has been given by Z. X. Yang. The cyclic representation was given 
in@ 

Heisenberg automorphism of 5 T r? 

By shifting rj — > rj + \[ol\t + a.2) and z\ — > z\ — h( a l T + a 2), from lemma 1 and lemma 
2 it is easy to check that there is an automorphism: 

S p —> S'p = S p+a uo a ^ 2+ ^ , F 7 —> F; = u a ^F 7 (73) 

which satisfy 

W = TiVpUK (74) 

b 

Remark: The automorphism for commutation rule of S a and similar change in basis for 



general level I representations have been given in the paper [ 40 1 
A n and su n (T) vs. S T)rj , su n (T) ~ lim^o $t, v 



Let rj — > 0, the S a (|62j), ( p3|) in 5 Tj ^ degenerates into the E a ( |32j ) and (|33D in su n (T) in 



the first order of 77, 5 a — > E a , (a / 0, 0) while So,o — ► Eqq is the zeros order. 
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The su n {T) is an ordinary Lie algebra with highest weight representation. But defining 
relations of S Tr) are quadratic, no linear sum of S a lies on the center of <S T „ and S Tr) has no 
highest weight representation. 

The S a is gauge transformed by the Z n x Z n Gji{n) shifts of Zi among coverings, and 
meantime itself transforms nonabelliantely from S a to S a +p by the G-^(n) shift rj — ► 77 + 
h{P\ + (hr)- The S a implicitely intertwin with the Z n x Z n symmetries. For example, in 
(I a/3 )ij, i,j are the indices of the solitons, a, (3 are the indices of branes (coverings). Thus the 
affine Weyl reflections consist of the reflections of solitons (in the center mass system) and 
the truncation by a period of T n into another covering. So the operators S ai0l2 interchange 
the wave functions from brane a± to a 2 . The change of the strenth B of the fields; i.e. the 



strenth r\ by „ induces the su(n) gauge transformation of the operators S a to S a +p. 
The Z n Z n Belavin model 
Let 



(75) 



where W a (u) = CT " ( " +r?) . Then 



(v) 

Ri 2 {u\ - U2)L 1 (u 1 )L 2 (u 2 ) = L 2 {u 2 )L 1 (u 1 )R^ 2 (u 1 - u 2 ), (76) 

where R B {u)f k = J2a W a (I a )l '(I^ is the Belavin Z n x Z n vertex model. Let rj — > 0, then 
R B — > R G . Furthermore following |41] using the same gauge transformation as we did from 



L G ^ L \rmCM 

one can obtain the Ruijsenaars operator TL n as the u independent part, this 
is very important to understand the dynamics of soliton. And the det \L B (u) — k\ will give 
the Ruijsenaars 7i n at the n-th power of k: det|L B (u) — k\ = gives the spectral curve. 
All this will be related to the elliptic quantum group E T ^ and to be discussed in the next 
section. And its relation to the N.C. field theory such as the N.C. CS and gauged WZW will 
be investigated in the future. 



6 Elliptic quantum group and IRF model 

To define the elliptic quantum group E TtTI (sl 2 ), we start from a n x n matrix R(u, A) G 
End(V (8> V), where V is the element in representation space of sl(n) and A a weight vector 
A = (Ai,A 2 ,---,A n ) 6 h*. 

n 

R(U, X)=J2 E hi ® E V + iZ a ( U ' X ij) E i,i ® E hi + J2 X ij) E i,3 ® E 3^ ( 7T ) 

where Ei^e^ = Sj,k&i, e k = (0, • • • , 1, • • • , 0), here 1 is at the k-th place, and 

e(u)e(X + r ] ) 9( U + \)e(rj) 

A) = 9(u- V )9(\y ft"' A) = 6{u - rj)0{\) ' 
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It satisfies the dynamical YBE: 

R(ui,u 2 , A - rjh®) 12 R(ui, X) 13 R{u 2 , A - r/h^) 23 

= R{u 2 , X) 23 R{ Ul , A - r]h (2) ) n R{ Ul - u 2 , A) 12 (78) 

where R(u, A — r/h^) 12 acts on a tensor v\ tg) v 2 (g> V3 as i?(u, A — 77/i) (8) Id if has weight /i. 

The elliptic Quantum group E T)r) {sl n ) is an algebra generated by a meromorphic function 
of a variable h and the matrix elements of a matrix L(z, A) G End{V) with noncommutative 
entries, subject to the following relation 

R(ut -u 2 ,\ — rjh^) 12 L( Ul , X) 13 L{u 2 , A - V h^) 23 

= L{u 2 , X) 23 L(ui,X - V h^) 13 R( Ul - u 2 , A) 12 . (79) 



On the other hand, there is the well-known interaction-round- a- face (IRF) model. We 

denote by W(a, 6, c, d, u) its Boltzmann weight for a state configuration ^ ^ J round a 

face with a,b,c,d G /i* which is the weight space of the Lie algebra sl n . The only nonvanishing 
Boltzmann weights [42] can be represented as 

9{u + 7]') 



W(m, m — ej,m — 2ej ,m — ej,u) 



9(u + r]'m jk ) 



W(m, m - ej, m-ej- e k , m - e j: u) = J — , (j / k) (80) 

9(r]'m jk ) 

W(m, m — ej,m — e ? - — e k , m — e k , u) 



9{rf) 
9(u)9(r]'m jk - jf) 



0{rf)9{im jk ) 

where rf = —r], rrij k = mj — m k . They satisfy the YBE: 

^2 W(a,b,c,g,ui - u 2 )W{g,c,d,e,ui)W{a,g,e, f > 2 ) 
9 

= W(b, c, d, g, u 2 )W(a, b, g, f, ut)W(f, g, d, e, m - u 2 ). (81) 
g 

The -R-operator as n x n matrix in auxiliary spaces becomes the Boltzmann weight of the 
IRF model W(a, b, c, d, z): 

— R(u, 7]m)ej (£> e k = W(m, m — ej,m — ej — e k ,m — e^, u)ej (£> ej, (82) 

where e, + = ej + e k . Thus, eq. (|79|) can be rewritten as 

£ W(a, b, c, d, ui - u 2 )L« (ui , A) J fe L( 2 ) («2, A - kr,%, 

■I 

]T (txa, A)JlW( ttl ,A - iir?)^(a', 6', c', cT, «i - u 2 ), (83) 



3,3 



11,31 
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where 

rja = X — T)fv® = A — rjfi, rja = A, 7/ = —rj, (84) 

a — b = ej, b — c = Bji, a — d = ep, d — c = ei, (85) 

a — b = efc, b' — c = e&/, a — d' = ej 1 , d — c = e^, (86) 

e i + e i' = e z + e Z'> e fc + e fc' = e n + e ii • ( 87 ) 

We will drop the superscripts 3 of the quantum spaces later. 

Compare (|8~T| ) and (83), using explicite (80), it is easy to give an evaluation representation 
of the elliptic quantum group 

y - _ <To(u + |- ^-7ya fc j-gf^) n (T (-| + W 



which can be factorized as 43] 



% 

-ejt' 



L( i )(, ) A)=^^(c ej ^+^: 

i 

where 

I 0, otherwise 

and 5 = X^o 1 a «' a = 7 7 X (A - t}h). 

The evaluation Verma module V\(u) of -EV n (sl n ) is acted by the L(tt, \) 3 k with £ = A, and 
a /(/»), 

f(h)L{ = L{f(h-U-k)). 

The basis of V\(u) satisfies 

f(h)e k = f(A-k)e k , (fc>0); 4e = 0, (j < fc). 

If A + 1 = n + a+ P T , then the vector eo is the highest one in V\, the subspace of V\(u) 
spanned by e^, (k > n — 1) is a submodule. The quotient space L\(u) is a module of dim n. 

Following p]J, for an E T ^(sl n ) module V, let -F un (V) be the space of meromorphic func- 
tions of A with values in V and 1 periodic, F(X + 1) = F(X). Introduce indomorphism f(h), 
L{(w) G End(F un (V)) by the rule 

(f(h)F)(X) = f(h)F(X), (90) 
and the difference matrix operators 

(Li(w)F)(X) = L{(w)F(X - (j - k) V ). (91) 
This defines a difference matrix operator algebra -A(V) on V. 
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Factorizations of vertex model and IRF model 

Substituting the S a (z) 

as £ a + 



[U, z 



into the Belavin's L B {u)\ fl75|), we find ||] it can be factorized 
~. Thus, the intertwinner in (|89|) 0„ ~ 0A-e intertwines the 



Heisenberg shift of — % in 



n 

1 
2 



for the i-th brane and the change of r\\ a (r]z a ) by r] 



times unit weight vector (z a — ^ Ylk z k)v f° r the a-th soliton. 

The Heisenberg shift corresponds to the W a matrix for wrap on torus T, while the later 
correspond to the difference of A (~ z) in sections of bundle C 

Transfer matrix and Ruijsenaars operators 

The transfer matrix is a difference operator acting on the space .F(V[0]) of meromorphic 
functions of A G h* with values in zero weight space V[0] of the quantum space (h module) 
V. It is defined by the formula 



T(u)f(X) = J2L u (u,X)f(X-r ] h). 



(92) 



On the otherhand, as shown by Hasegawa by expanding the normal ordering det \ L B (u) 
1 41], one obtain the Ruijsenaars Macdonald operator M 



M 



N 

En 

% j-.j^i 



9{\j - Xj + Irj) 
0(X l - Xj) 



Ti, 



(93) 



and the other higher Hamiltonian. 

The Ruijsenaars operators (93) are integrable difference operators that give a q- deforma- 
tion of the Calogero-Moser integrable differential operators (|58|), inversely Calogero-Moser is 
the 7] — > limit of Ruijsenaars- Schneider. 



Let us as [44] consider the transfer matrix associated to the symmetric power S m V . The 
zero weight space of this module is trivial unless m is a multiple of n. If m = nl then the 
zero weight space is one-dimensional and is spanned by the sum of the tensors (g> ■ ■ ■ (8) e. 



over all sequences (ij) such that each number between 1 and n occurs precisely I times. Let 
us denote this sum by e. 



Let us as [45| identify S nl (C N )[0] using the basis e, and let T(z) be the transfer matrix 



associated to S n V(0). Then 



T(u) 



9{u - rjl) 



M 



(94) 



9{u — rjnl) 

which is the same as from f43fl by choosing the parameters 5, £ appropriately. 

Let us consider the I = 1 case, the Aj ~ Zi is the position of our solitons. All the M 
and other higher Hamiltonian depends only on Zij, and independent of neither spectral u nor 
the center Yl z i- So the long distance interactions, the dynamics of solitons are given by the 
elliptic Ruijsenaars-Schneider models. 
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Bethe Ansatz for the Hilbert space of non-commutative torus 

Now we turn to find the common eigenvectors of the transfer matrix and the Macdnoald 
operator M (p3|). It is easy to check that the vector v a (^) in finite dimensional space V n are 
eigenvectors of M (3!|. But the vector space V n as the trivial U(n) boudle |l5|] transformed 
only by I a P matrix, representing the Wi, independent of the twist i]. 

The noncommutative algebra A^{n) for the twisted case should be acted on the Hilbert 
space Ti„(n) ~ C<g>V n constituted by vector- valued functions rather than on finite dimensional 



vector space V n . The contribution of paper [21| on the algebraic Bethe ansatz for the elliptic 



quantum group is just to use transfer matrices acted on spaces TL n of vector-valued functions. 

For the module 0] L^(u) (= H n ) (we restricted this irre. case), the transfer matrices 
T(u) preserves the space F un (H)[0] of functions with values in the zero- weight space H[0] 
and commute pair wise on it. The highest weight vector of H is a v £ F un {7i)(h) such that 
Lj(u)v = (j > i), and F un (TC) is spanned by the vector of the form 

Lf 1 (h)---L^(t m )v (95) 

with all j m > i m . 

Its eigenfunction is similar to that for the elliptic Gaudin. Now the Hermite Bethe pa- 



rameters satisfy the so called nested Bethe ansatz equation |47|, 



The Bethe ansatz equation (BAE) (|96|) is the condition of t\,---,t m and v = g(X)eo, 
so that (|95| ) is an eigenvector. The BAE is determined by the vanishing of the unwanted 
terms [^]. When the diagonal L^,(w) in T(w) move through the lowering L^(t) to v, the 
contribution of different diagonal L should be cancelled. 

The noncommutativity causes the shift of Aj by rj. So besides the usual BA. condition 
for ti and Aj, we following [^] (cf. their proof of Lemma 3 for detail), find that as the / = 1 
case 

Check carefully for the vanishing of unwanted terms, the second factor of g(\) ensures 
that the first A dependent factor of the V n matrices L\ and L 3 - (|8|) scattered through L*- 
be cancelled. Meanwhile the difference action on C{u) of L\ and shift the e aX factor 
differently, with a phase shift , which is independent of A. This independence of r] is 
sufficient to ensure the usual BA. eq. for different ij. But meanwhile the wave function will 



have the form [21 



v;=n eCiAi ip( A *+^-'?)> ( 98 ) 

i j=l 



22 



which will be twisted by Qr] when Aj changed by the difference action of noncommutative 
Wilson loop U\, U2 i.e. by diagonal L\ and nondiagonal L? k respectively. 

In case of XYZ or equivalently the sl(2) IRF, in the Gaudin limit [28], the parameter 



(j = 2nvj, (0 < nuj < is related to the twist in the boundary actions [49]. It describes 
the twist of the wave functions |)0], |5l]] 

^ = e 2 ^ v+2av . (99) 

aeZ 

To ensure the convergence of the series in case r\ equals some rational time of 1+ ™ T , choose 
v = — , then the Fourier series truncated to finite n terms [50], and we will obtain cyclic 
representation as describe in pq]. 

For the r\ general case, the zero- weight condition restricted the number m of Lj*(ti), which 
will be the same as zero-weight space nl, (in our case I = 1). In case n = 2, this implies 
the number of lattice N should be 2 (n = 2) times the number of operators Lj (their Bj). 
The above are arguments in usual algebraic B.A. on vector space V n without the difference 
operator on C(u), which we turn now. Now the finiteness of the dimension is the truncation 



of the Verma module to its quotient space [ 46 1 



As our discussion for QHE, we choose Q to be an integer A, (we have normarlized A to 
be 1, i.e. case I = 1), this will give us the wanted noncommutativity. Actually this is just 
the contribution of difference operatoin when TJ\ like L|, L^t^ commute with U2 like L\ +l . 

7 Discussion 

There is a lot to be investigated later, e.g.: How to use the Weyl Moyal transformation 
to give the explicit functional expressions of the soliton solution i.e., the projection and ABS 
operators; how the affine Weyl reflections imply for the duality and Morita equivalence. It is 
easy to find ADE orbifold and BC orientfold case for Gaudin and CM. models. 

For the rj ^ case, the N.C. q KZB and q Virasoro should be considered. The N.C. R.S. 
as C.M. can be obtained as Hitchin systems from the cotagent bundles on N.C. torus. 

Thus in anticipation for the application in N.C. models, we collecte some known elliptic 
algebra in the second part of this paper. 

The Sklyanin algebra is the deformed su n (T). (The differencial operator is deformed into 
difference of rj ) . The Lax pair and i?-matrix of Z n (g> Z n Belavin is the deformed L-matrix and 
r-matrix of elliptic Gaudin. The cotangent bundle on torus for algebra will be deformed into 
that for group |5^ ]. The C.M will be deformed into RS related with IRF and the quantum 
group. 



The corresponding field are the gauge WZW and CS [53], the N.C. case has been studied 
recently by [54]. 



We will investigate more carefully N.C. aspect of the integrable many-body and gauge 
field [55] and string brane [^] later. The dynamical variables 2«(~ Aj) and <9j in exponents 
will be the z% and Zj(~ di) for N.C. soliton. The N.C. source now equals the U1U2U1 as 
the affine Heisenberg Double [57] and the N.C. Wilson loop for the RS reduced form CS 54]. 
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